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Abstract 

The effects of nuclear re-interactions in quasi elastic and resonant neutrino 
interactions have been considered in the framework of the nuclear models of 
the DPMJET code. A preliminary investigation on the modifications induced 
on the final state has been performed. Some consequences affecting the ex- 
perimental identification are discussed. 
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1 Introduction 



The interest in neutrino interaction is nowadays mostly oriented to the question 
of neutrino mass and flavour oscillations. The anomaly in atmospheric neutrino 
flux, pointed out by many experiments [[I], @, |3|, suggested the possibility to explore 
the region Am 2 ~ 10 _2 -10 -3 eV 2 , using a Long Base Line (LBL) neutrino beam. 
Considering the recent CHOOZ result 0], ruling out the v^— > v e explanation for the 
atmospheric anomaly, the interest for LBL experiments is shifted to v T appearance 
search with energy £^4-10 GeV, since the recent result from Superkamiokandeg], 
folded with the older Kamiokande and 1MB results, suggests oscillation parameters 
sm 2 2#~l and Am 2 ~ 5-10~ 3 . The v T tagging is rather difficult at these energies 
both from a statistical and an experimental point of view. The first point comes 
from the overlap of the low cross section at these energies with the poor neutrino 
flux in the LBL case, with respect to a Short Base Line situation, decreasing like ~ 
r -2 , where r is the distance between the neutrino production point and the neutrino 
detection point. The second point is related to the intrinsic difficulties in measuring 
the kinematic variables of v T interaction products at low energy. 

However, at these neutrino energies, the selection of quasi elastic and resonant 
events can be useful to identify v T candidates. At energies not far from the thresh- 
old for the production of r lepton, the quasi elastic cross section gives still a large 
contribution to the the total charged current cross section. Moreover the intrinsic 
simplicity of the kinematics, simplifies the search "a la Nomad" || for v T identifica- 
tion. Favourable are those events in which the r lepton decays into an electron or 
muon (plus neutrinos) or into a single pion (plus neutrino). 

However, the intrinsic cleanness offered by kinematics of the quasi-elastic inter- 
action of neutrinos on single nucleons, can be obscured in case of nuclear targets. 
Nuclear effects, like Fermi motion and nuclear re-interaction of the nucleon inside 
the nucleus, have to be considered with care. These effects can provide momentum 
imbalance, the production of additional particles and may lead to the misidenti- 
fication of the outgoing nucleon. All this phenomenology can deteriorate the sig- 
nal/background ratio, the event recognition efficiency and the vertex identification. 
Moreover, in case of experiments that plan to use emulsions as active targets, one 
expects to see, sometimes, low energy particles and fragments close to the interaction 
vertex ( "grey" and "black" particles, in the language of emulsions) . 

Guided by this consideration we start from the quasi-elastic and resonant neu- 
trino scattering simulation code of [|7|], which already takes into account some nu- 
clear effects, like Fermi motion and Pauli blocking. We have interfaced it to the 
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DPMJET-II-.4.1 code || ^, |T0], which was created to treat hadron-nucleus and 
nucleus-nucleus interaction at high energy but which also contains a Formation 
Zone Intranuclear Cascade (FZIC) for low energy interaction of the produced sec- 
ondaries with the spectator nucleons, nuclear evaporation and the formation of the 
residual nucleus. The nuclear models of DPMJET have been tested against a wide 
set of experimental data. In particular, the main check of the de-excitation and 



fragmentation algorithms comes from the comparison with emulsion data [11, 12 



More refined models which take into account the effect of the nuclear poten- 
tial on the particle trajectories exist 0, pj], [15|, 16, 17 and will be considered by 



other authors |18| also in the framework of neutrino simulation. While waiting for 
these more refined codes, we think that we can obtain from the use of DPMJET a 
lot of essential information, which, as a first approximation, is reasonable and not 
achievable from presently existing codes. 

In this paper we present the first results obtained using this code, on the com- 
parison with single nucleon interaction with and without considering the nuclear 
effects and the implication for v T explicit search. 



2 Description of the code 

The code extends the quasi-elastic neutrino-nucleon and neutrino-nucleus model 
QEL0, based on the formulation of Lewellyn-Smith[ |19f| , to neutrino collisions on 
nuclear targets including the intranuclear cascade and nuclear fragmentation and 
evaporation. The QEL code has been transformed to double precision, to match the 
DPMJET environment, and it has been modified to use the Fermi momenta of the 
nucleons from DPMJET. 

All neutrino flavours can be considered. In the case of tau lepton production, the 
leptonic decay into muon (or electron) plus neutrinos can be optionally activated. 
For this reason, we have also introduced the calculation of the polarization of the 
final lepton, according to the treatment of EO . 



The nucleon generated by the QEL code can re-interact in DPMJET inside the 
nucleus according to the Formation Zone Intranuclear Cascade (FZIC) model [EI], 



22] , |Tlf contained inside DPMJET. Secondaries from this first collision are followed 
along straight trajectories and may induce in turn intranuclear cascade processes if 
they reach the end of their formation zone inside the target, otherwise they leave 
the nucleus without interaction. Inelastic secondary interactions in the FZIC are 
described by the Monte Carlo code HADRIN p3l |2^|. For the sampling of elastic 
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nucleon-nucleon scattering below 4 GeV the parameterization of the HETC-KFA 



code |26| was adopted. 

The treatment of nuclear effects within the MC model has already been discussed 
in [27] and in more detail in [11, 12 1. Since they are essential in calculating excitation 
energies of nuclei left after primary interactions and intranuclear cascade processes 
we summarize the basic ideas. Fermi momenta for nucleons as well as a simplified 
treatment of the nuclear potential are applied to control the generation of low-energy 
particles. Nucleon momenta are sampled from zero-temperature Fermi distributions 

a = N n ' p P (1) 
dp (p n F ' p ) 3 { ' 

The maximum allowed Fermi momenta of neutrons and protons are 



n,p 

Pf 



iV n ' p \ 3h 3 



V A J 8vr 



(2) 



with V A being the volume of the corresponding nucleus with an approximate nuclear 
radius Ra = r ^4 1//3 ,r = 1.29 fm. 

Modifications of the actual nucleon momentum distribution, as they would arise, 
for instance, taking the reduced density and momenta in the nuclear skin into con- 
sideration, effectively result in a reduction of the Fermi momenta as compared to 
those sampled from Eq. (H). This effect can be estimated by a correction factor a^ od 
which modifies the Fermi-momenta. Results presented in this paper have been ob- 
tained with a^ od =0.60. The depth of the nuclear potential is assumed to be the 
Fermi energy and the binding energy for outer shell nucleons 

V»> p = + EZ d - (3) 

To extend the applicability of the model to the energy region well below 1 GeV an 
approximate treatment of the Coulomb-potential is provided. The Coulomb-barrier 
modifying the nuclear potential is calculated from 

Vc = -* ZA , (4) 

with the mass numbers A\,A2 and charges Z\^Zi of the colliding nuclei, i.e. with 
A\ — \Z\\ — 1 for charged hadrons entering or leaving the target nucleus, e denotes 
the elementary charge and r = 1.29 fm. 
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The excitation energy U of the residual nucleus with mass number A Tes and 
charge Z ies , i.e. the energy above the ground state mass E 0tres , is given as 

U = E -En 

w -'-'res - LJ 0,resi 

Eq,tw = ■ZresW'p + (Ares — Z ves )m n — i?bind (Aes > ^res)- (5) 

The binding energy E hind (A res , Z res ) is obtained using the experimentally determined 
excess masses of all known (measured) nuclides and using mass formulae for nuclides 
far from the stable region, where no measurements are available. The excitation 
energy is obtained within our model from an explicit consideration of the effects 
of the nuclear potential (Eq. ([$])) and the Coulomb energy (Eq. (f|)), i.e. from 
corrections which are applied to the 4-momenta of the final state hadrons leaving the 
spectator nucleus. We modify the energies of these hadrons by the potential barrier 
and rescale the 3-momenta correspondingly. It is assumed that these corrections 
have to be applied to nucleons wounded in primary and secondary interactions and 
only to those hadrons which are formed inside the spectator nucleus corresponding to 
the sampled formation time. Among these particles we find apart from the nucleons 
a small fraction of other baryons, which are assumed to move in a nucleon potential 
and mesons to which we apply an effective meson potential of 0.002 GeV. Due to 
energy-momentum conservation these corrections lead to a recoil momentum and, 
therefore, to an excitation of the residual nucleus. In addition, there is a further 
contribution to the recoil momentum of the residual nucleus arising from potential 
corrections applied to the momentum of the projectile hadron entering the nuclear 
potential and from cascade nucleons with kinetic energies below the nuclear potential 
which are therefore not able to escape the spectator nucleus. Pauli blocking is 
considered, and events are rejected accordingly. Instead, no angular momentum 
barriers are considered. 

At the end of the intranuclear cascade the residual nucleus is supposed to be 
left in an equilibrium state, in which the excitation energy U is shared by a large 
number of nucleons. Such an equilibrated prefragment nucleus is supposed to be 
characterized by its mass, charge, and excitation energy with no further memory 
of the steps which led to its formation. The excitation energy can be higher than 
the separation energy, thus nucleons and light fragments (a,d, 3 H, 3 He) can still be 
emitted: they constitute the low-energy (and most abundant) part of the emitted 
particles in the rest system of the residual nucleus, having an average energy of few 
MeV. 

In heavy nuclei the evaporative process is in competition with another equilib- 
rium process, that is fission p5[ . For the fission probability, a statistical method 
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can be used [^, [I(J . 

Other de-excitation mechanisms are more suitable for light residual nuclei. The 
one adopted for this calculations is the so called Fermi Break-up model |JT], [52 



where the excited nucleus is supposed to disassemble just in one step into two or 
more fragments, with branching given by plain phase space considerations. Accord- 
ing to the picture of the compound nucleus like an equilibrated system determined 
only by its mass, charge and excitation energy, with no memory of previous steps 
of the interaction, Fermi Break-up is activated in the model every time the current 
compound nucleus has mass number A < 17, including possible light fission frag- 
ments. The fragmentation of higher mass compound nuclei is not yet included in 
the model. This process, although its cross section is quite small, is important when 
considering the distribution of residual nuclei, because it can produce isotopes very 
far both from the target mass and from the fission product distribution. 

The evaporation stage ends when the nuclear excitation energy becomes lower 
than all separation energies for nucleons and fragments. This residual excitation 
energy is then dissipated through emission of photons. 7-de-excitation proceeds 
through a cascade of consecutive photon emissions, until the ground state is reached. 
The cascade is assumed to be statistical as long as the excitation energy is high 
enough to allow the definition of a continuous nuclear level density. Below a (some- 
what arbitrary) threshold, set at the pairing gap value, the cascade goes through 
transitions between discrete levels. In reality, photon emission occurs even during 
the preequilibrium and evaporation stages, in competition with particle emission, 
but its relative probability is low, and it is presently neglected in the model. 



3 Simulation results 

In order to test the features of the code, we considered 5,000 quasi elastic interactions 
of 10 GeV on light (Carbon), intermediate (Silicon) and heavy (Iron) targets. We 
immediately see that additional particles appear in the final state. Table 1 shows 
the comparison of the results on different targets, as far as the average nucleon, 
gamma and charged pion multiplicities and the average de-excitation gamma en- 
ergy are concerned. Moreover, we show the probability P to obtain a "clean" two 
body final state (just muon plus proton) and the average missing P±. We notice 
that the main effect of the nuclear re-interaction consists in additional nucleon pro- 
duction, although at very low momentum, as it will be shown later. Charged pion 
production(whose kinetic energy is around 300 MeV) is rather negligible. 
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Target 


< N n , p > 


< N w ± > 


< N^> 


< E 1 > 


P(n=2) 


^ pmiss ^ 






(MeV) 




(MeV/c) 


free nucleons 


1. 











100% 





Carbon 


1.49 ±0.02 


0.016 ±0.002 


0.79 ±0.02 


2.1 


27.4% 


175 


Silicon 


2.02 ±0.02 


0.025 ±0.002 


2.68 ±0.02 


1.5 


0.% 


183 


Iron 


2.26 ±0.05 


0.023 ±0.005 


3.42 ±0.06 


1.5 


0. % 


216 



Table 1: Comparison of the result of interaction (E y ^=10 GeV) on different 
targets. 
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Figure 1: Example of a simulated quasi elastic event at 10 GeV in a glass- 
scintillator sampling calorimeter (nuclear re-interaction off). 
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To give an idea of the event topology distortion introduced by nuclear effects, we 
simulated a fine grain low density sampling calorimeter (1/4 Xo glass + 1.5 cm of 
liquid scintillator) using the FLUKA package^]. Fig. 1 and 2 show the same QEL 
event, without nuclear effects (Fig. 1) and and with nuclear effects (Fig. 2). It's 
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Figure 2: The same event of Fig. [l| when nuclear re-interactions are considered. 

evident, looking at these figures, that the identification and tracking of the proton 
becomes more difficult, both for the presence of extra particle in the final state and 
for the modification of the kinematics of the proton itself. 

Fig. 3 shows the momentum of the protons produced in the interaction. In 
the top part we consider all the protons emitted while in the bottom part, we 
considered only the event leading proton (the most energetic one). We show in the 
same figure the case with and without nuclear re-interaction. We learn that the bulk 
of the additional protons is produced with low momentum, as typical of intranuclear 
cascade products. As expected, the effect tends to become more significant for heavy 
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Figure 3: Protons (top) and leading proton(bottom) momentum distribution, in 
f M +n, f M +C and f M +Fe interaction (E u = 10 GeV) 

target nuclei. There is also some dependence on energy. The effects of nuclear re- 
interaction manifest themselves in the distortion of the kinematic variables. As an 
example, in Fig. 4 we show the distribution of the angle between the muon and 
the outgoing leading proton. The case with and without nuclear re-interactions are 
distinguished and compared to the case of the free neutron target. We notice how 
the proton now, even the leading one, can access much larger angle with respect 
to the simple kinematics of pure neutrino-neutron quasi-elastic scattering. Further 
kinematic effect can be observed looking at event missing P±. Also in naive models 
we expect the appearance of non zero missing P± from Fermi motion. Nuclear re- 
interactions increase this effect ( Fig. 5 and 6 for C and Fe target respectively). 
There, we show in the top part the change in missing P± distribution from the case 
without intranuclear cascade (left) to the case including re-interactions (right). A 
tail in the last case is evident extending at large value of P±; this is a dangerous 
effect for the experiments aiming to tag the v T using kinematic cuts "a la Nomad" || . 
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Figure 4: Distribution of the angle between the outgoing proton and the neutrino 
direction, in i^+Fe interaction, E v = 10 GeV) 

This is clear in the bottom part of Fig. 5 and 6, where we show the scatter plot 
of the angle between the proton transverse momentum P\_ and the muon P±(4>2), 
versus the angle between P± and the missing transverse momentum p™ lss (<pi). For 
completeness we show in Fig. 7 the case of z/ T +n^r+p, where we considered the 
channel r— tfi+v+i? in Fe target. Due to nuclear effects, we observe the presence of 
events in which interactions give small </>i angles, close to the region expected for 
v T candidates. 

The intranuclear cascade is relevant as far as the and v e are concerned, while 
in this last case (v t ,t—*ij,+v+P) it is obscured by the missing P± introduced by the 
t decay. In any case, it seems that, despite the smearing introduced by nuclear re- 
interactions, a region still exists in the 4>2-(j>i plane which allows the identification 
of r candidates. Just to consider an example, although only at particle level, if we 
take , the (non optimized) region on the right of the triangle drawn in the rightmost 
bottom plot of Fig. 7, we find that the efficiency of detecting a r candidate is 
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Figure 5: z^+C interaction^^ =10 GeV). Missing P± distribution (top) without 
(left) and with nuclear re-interaction (right). Scatter plot of <pi and cf>2 angles 
(bottom) without (left) and with (right) nuclear re-interaction. 

864/4584, having rejected 100% of 4514 interactions. Of course, this is just an 
over simplification, since in practice all the other resolution effects coming from a 
realistic experimental simulation have to be considered. 

4 Delta resonance 

Delta resonance process has been included in our code, for both CC and NC interac- 
tions. As a relevant example for experimental application, we have considered 5,000 
CC interations of 10 GeV on Iron. All the channels A ++ ^p+7r + , A + ^p+7r° and 
A+^n+7r + have been considered, and the results are shown in Table 2. As in the 
case of QEL interaction, we see the presence of extra particles in the final state: in 
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Figure 6: v^+Fe interaction^^ =10 GeV. Missing P± distribution(top) without 
(left) and with nuclear re-interaction (right). Scatter plot of <pi and 4>2 an- 
gles(bottom) without (left) and with (right) nuclear re-interaction. 

practice, for Fe target, we never observe events preserving the final state obtainable 
without nuclear reinteractions. It might be intersting to notice how in a non neg- 
ligible fraction of cases the charged pion can be absorbed inside the nucleus. This 
complicates the separation of quasi elastic interactions from resonance excitation: 
these two classes must be always considered together. Also, the effect of the nuclear 
rescattering manifests in an additional tail of event missing momentum. 



5 Conclusions 

We have shown how the nuclear re-interaction in quasi elastic neutrino events can 
play an important role in the experimental neutrino detection. The nuclear models 
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Figure 7: is T +C interaction^;, =10 GeV. Missing P± distribution(top) without 
(left) and with nuclear reinteraction(right). Scatter plot of (f>\ and 02 angles(bottom) 
without (left) and with (right) nuclear reinteraction. A possible selection area to 
identify r candidates is the one to the right of the triangle drawn in the last plot. 



contained in the DPMJET code seems to be adequate for the description of the 
phenomenology of neutrino interaction in view of design of future detectors for Long 
Base Line experiments. Of corse this is just a preliminary stage of the work. Further 
improvements will allow the inclusion of scattering. Nuclear effect introduces also 



changes in the total neutrino-nucleon cross section, as discussed in [33 
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Process 


<N P > 


< N n > 


< N n ± > 


< N n o > 


< N^> 




2.14 


1.97 


0.85 


0.03 


3.68 




2.13 


1.57 


0.05 


1.02 


3.46 


A + ^ri7r + 


1.64 


2.65 


0.91 


0.04 


3.39 



Table 2: Comparison of the result of interaction ( E v =10 GeV ) on Iron for 
different resonant processes 
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